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We present a Monte Carlo renormalisation group study of the SU(2) gauge theory with two Dirac 
fermions in the adjoint representation. Using the two lattice matching technique recently advo- 
cated and exploited in Ref. [1], we measure the running of the coupling and the anomalous mass 
dimension. 
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1. Introduction 

Technicolor theories with fermions in higher representations of the gauge group can poten- 
tially provide a dynamical electroweak symmetry breaking mechanism without conflicting with 
electroweak precision data. Minimal Walking Technicolor is an example of such a theory, a SU(2) 
gauge theory with two Dirac fermions in the adjoint representation [2, 3], It is expected from 
perturbation theory to be in or near to the conformal window, although since this is at relatively 
strong coupling non-perturbative results are desirable. Lattice simulations have indeed found that 
the gauge coupling runs very slowly [4, 5, 6], more slowly than the perturbative prediction, al- 
though distinguishing between conformal and near-conformal behaviour is an inherently difficult 
task. In order to be phenomenologically viable the theory must have a large anomalous dimension 
(y ~ 1) [7, 8, 9]. The conjectured all-order beta function predicts y = 0.75 [10], and the anomalous 
dimension has been measured non-perturbatively in some recent lattice studies [4, 11], which have 
all found lower values. 

In this work we measure the anomalous mass dimension using the Monte Carlo Renormali- 
sation Group, a technique which has recently been used to investigate theories with many flavours 
of fermions in the fundamental representation [1, 12]. We consider the evolution of all possible 
couplings of the system under Wilson RG block transformations, where with each blocking step 
ultraviolet fluctuations are integrated out. Fixed points are characterised by the number of relevant 
couplings, which have positive scaling dimensions and flow away from the fixed point. Irrelevant 
couplings have negative scaling dimensions, and flow towards the fixed point, so that their IR values 
are independent of their UV values. 

2. Method 

With each RG step, changing the scale by a factor s, irrelevant couplings will flow towards 
the fixed point, and relevant couplings will flow away from it. After a few steps the irrelevant 
couplings should die out, leaving the flow following the unique renormalised trajectory (RT). If 
we can identify two sets of couplings which end up at the same point along the RT after the same 
number of steps, then they must have the same lattice correlation lengths, f = |'. Since the physical 
correlation length £, = qa should not be changed by the RG transform, this means that they both 
must have the same lattice spacing a, or inverse cutoff A -1 ~ a. If they end up at the same point, 
but one takes an extra step, then their lattice correlation lengths must differ by a factor s, and hence 
so must their UV cutoffs. 

To identify such a pair of couplings, we need to show that after n and (n — 1) RG steps re- 
spectively their actions are identical. Explicitly calculating the actions would be complicated, but 
instead the gauge configurations themselves can be RG block transformed, and showing that the 
expectation values of all observables on these gauge configurations agree is equivalent to directly 
comparing the actions that generated them. 

2.1 2-Lattice Matching Procedure 

Starting with the SU(2) pure gauge theory, where the gauge coupling is the only revelant 
parameter, the procedure is as follows. 
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1. Generate an ensemble of gauge configurations with an action 5(j8) on a L 4 lattice. 

2. Block these n times to produce an ensemble of configurations on a (L/s n ) 4 lattice, and mea- 
sure the expectation values of various observables on them. 

3. Generate a new ensemble of gauge configurations with an action S(jS') on a (L/s) lattice, 
for a range of values of /3'. 

4. Block each of these n — 1 times to produce an ensemble of configurations on a (L/s") 4 lattice, 
and measure the same observables for each j3'. 

5. Interpolate in j3' such that each observable after taking n steps on the larger lattice agrees 
with the same observable after taking (n — 1) steps on the smaller lattice. 

6. Repeat for different n, e.g. for s = 2,L = 32, three values can be used: n = 2,3,4. 

We have now identified, for each n, a pair of bare gauge couplings (j8,j8'), with lattice cor- 
relation lengths that differ by a factor s, | ' = In the limit — >• °°, the quantity A/3 = /3 — 
/3' = Sb{P;s) is the step scaling function for the bare gauge coupling. This is the analog of the 
Schrodinger Functional step scaling function for the renormalised coupling, a(u, s), and in the UV 
limit where g 2 — > = 2N / j3, there is a simple relation between the two: 

Sbjfijs) = a(u,s) _ 1 
j5 u 

We use the following s = 2 RG blocking transform [1] 



V n ,u = Proj 



(1 — Ot^Uf^nUfi+n^ + — U n yU n ^ v ^U n -\-^J r y^JJ n ^2^ l 



(2.2) 



Here a is a free parameter, which can be varied to optimise the transformation. Changing a 
changes the location of the fixed point, and how quickly the unique RT is approached in a given 
number of steps. Ideally it should be chosen such that 

• All operators predict the same matching coupling between (n,n — 1) pairs for a given block- 
ing step n. (Deviations are a measure of the systematic error from not being at exactly the 
same point along the RT) 

• Consecutive blocking steps predict the same matching coupling, i.e. the coupling for which 
(n,n— 1) pairs agree should be the same for all n. (Deviations show that the RT is still being 
approached in the irrelevant directions) 



3. Pure Gauge Results 

As an initial test of the method we determined the step scaling of the bare coupling for the 
SU(2) pure gauge theory. We matched in /3, using 200 SU(2) pure gauge configurations for each j8 
on 32 4 and 16 4 lattices. We matched in the plaquette and in the three six-link loops. An example 
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of the matching of the plaquette is shown in Figure 1(a). The red, green and blue horizontal 
lines show the average plaquette on the 32 4 lattice after 2, 3 and 4 blocking steps respectively, at 
/3 = 3.0, a = 0.57. The interpolated red, green and blue points show the average plaquette on the 
16 4 lattice after 1, 2 and 3 blocking steps respectively, as a function of /3', also at a = 0.57. The 
value of j8' where the two red lines intersect gives the matching coupling for n = 2, similarly the 
green and blue lines give the matching coupling for n = 3 and 4. 

This matching is repeated for each observable, and the spread of predicted matchings for each 
n gives a systematic error on the central matching value. The whole procedure is then repeated 
for various values of a, as shown in Figure 1(b). In the limit n — > oo the predicted matching value 
should converge on a single value that does not change if n is increased, and is independent of a. 
With a finite number of steps we can exploit the a-dependence to find an optimal value of a, where 
subsequent RG steps predict the same matching value. The intersection of the last two blocking 
steps gives a central value for 57, Q3 = 3.0; s = 2), while the range of couplings for which any of the 
blocking steps intersect within errors gives the uncertainty on this central value. This was repeated 
for other values of j8, and the resulting step scaling of the bare coupling is shown in Figure 2. It 
agrees well with the perturbative prediction. 



SU(2) Pure Gauge; beta=3.0; alpha=0.57; 32(1 6) plaquette matching beta-beta' vs alpha, 32(1 6) lattice matching, beta=3.0 




beta' alpha 

(a) Plaquette Matching (b) a-Optimisation 



Figure 1: An example of the matching of the plaquette in /3 for the pure gauge case. This is repeated for 
each observable to give a systematic error for each matching, then a is varied such that all blocking steps 
predict the same matching. 

4. Anomalous Dimension Results 

At an IRFP in the full theory with dynamical fermions, the gauge coupling is expected to be 
irrelevant, leaving the mass as the only relevant operator. This means we should in principle be 
able to match in the mass at arbitrary couplings, as long as we have sufficient RG steps for the 
gauge coupling to flow to its fixed point value. We match observables as for the pure gauge case, 
but instead of matching in /3, we set j8' = j8, and match pairs of bare masses (m,m'). 

We use the HiRep RHMC code [13] with a Wilson plaquette gauge action and adjoint Wilson 
fermions. We generated ~ 1000 configurations on 16 4 and 8 4 lattices, for a range of bare masses 
at each /3. This allows two matching steps, after 2(1) and 3(2) steps on the 16 4 (8 4 ) lattices. An 
example of this matching and subsequent a-optimisation is shown in Figure 3. 
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Figure 2: Bare step scaling Sf,(l5;s = 2) for the pure gauge theory. 



The bare mass is additively renormalised, so we convert the bare masses to PCAC masses. We 
measure the PCAC mass averaged over the two central time slices, in, as a function of bare mass 
for each j3,L, as shown in Figure 4, and use this to convert the pairs of bare masses to pairs of 
PCAC masses. 

— i r (l i\ , / i \ iH + ^ *)/a(^o) 

m = j [mpcAC ( 2 — 1 ) +m PC AC ml , m PC Ac{xo) = ttfi \ t 4 - 1 ) 

The anomalous mass dimension can then be extracted from the relation m' /rh = 2 r+l . If we 
are at an IRFP then this is a scheme-independent quantity, so we should be able to repeat the 
measurement at any value of j3 and find the same anomalous dimension. We used three values 
of j3, jS = 2.25,2.50,3.00, and all give consistent predictions for the anomalous dimension. The 
matching PCAC mass pairs and a linear fit are shown in Figure 5, and predict y = 0.49(13). 



5. Coupling 

It would also be interesting to measure the step scaling of the bare coupling s/, (jS ; s) for the 
full theory with dynamical fermions, as was done for the pure gauge case. A fixed point would be 
indicated by a change of sign in this quantity as the bare coupling is varied from weak to strong 
coupling. The difficulty is that the mass is a relevant operator, while the coupling is expected to be 
at best nearly marginal, so that in order for the MCRG to pick out the behaviour of the coupling 
the mass would have to be tuned to zero. Furthermore, even if the mass is tuned sufficiently close 
to zero that we are initially following the evolution of the gauge coupling (assuming it is the least 
irrelevant remaining operator), we can no longer take the limit n — > °°, because the coupling will 
flow to its fixed point value, and the flow of the mass will eventually dominate. 
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SU(2) Adjoint; beta=2.50; alpha=0.65; m=-1.05; 16(8) plaquette matching 

38 « 



Ipha. 16(8) lattice matching, beta=2.50,m=-1 .05 





0.72 0.74 



(a) Plaquette Matching 



(b) a-Optimisation 



Figure 3: An example of the matching of the plaquette in the bare mass m. This is repeated for each 
observable to give a systematic error for each matching, then a is varied such that all blocking steps predict 
the same matching. 



PCAC mass L=8 PCAC mass L=1 6 




(a) L = 8 (b) L = 16 



Figure 4: PCAC mass as a function of the bare mass on 8 4 and 16 4 lattices for J3 = 2.25,2.50,3.00 



The method would therefore be to tune the mass close to zero initially, then take a few RG steps 
where the flow is following the gauge coupling, and extract the running of the coupling before the 
flow in the mass becomes significant. This has been successfully performed in Ref. [12] using 
staggered fermions for SU(3) gauge theories with 8 and 12 flavours of fundamental fermions. 

6. Conclusion 

We measure the anomalous mass dimension for Minimal Walking Technicolor and find y = 
0.49(13), in line with previous lattice determinations, and smaller than would be desired for phe- 
nomenology. A large anomalous dimension y = 1 is strongly excluded by our data, with a % 2 jd.o.f 
of 8.7 with 12 degrees of freedom. 

We intend to repeat this calculation on larger lattices to reduce our systematic errors, and with 
more observables to improve our determination of the systematic errors themselves, as well as 
with different RG blocking transforms to check for universality, and to allow us to go to stronger 
coupling. We are also planning to match in /3 with all bare masses tuned such that the PCAC 
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Figure 5: Matching mass pairs. Linear fit gives 7=0.49(13). Also shown is the linear fit given by assuming 
7=1, such a large anomalous dimension is clearly excluded by our data. 

masses are as close as possible to zero, in order to look for a fixed point by directly measuring the 
running of the coupling. 
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